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Chiral Magnetic Effect from Q-balls 
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We apply a generic framework of linear sigma models for revealing a mechanism of the mysterious 
phenomenon, the chiral magnetic effect, in quark-gluon plasma. An electric current arises along a 
background magnetic field, which is given rise to by Q-balls (non-topological solitons) of the linear 
sigma model with axial anomaly. We find additional alternating current due to quark mass terms. 
The hadronic Q-balls, baby boson stars, may be created in heavy-ion collisions. 
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It is widely believed that QCD has a phase transition 
between the hadronic phase and the quark-gluon plasma 
(QGP) phase at finite temperature and density. Exper- 
imental searches for QGP in relativistic heavy ion colli- 
sions have been revealed that QGP has highly nontrivial 
properties, such as its perfect fluidity [lJ4a, for exam- 
ple. The chiral magnetic effect (CME) [J, J5| is one of 
the most striking phenomena in QGP which has been 
recently studied from the theoretical and experimental 
viewpoints. The CME, the separation of electric charge 
along the axis of an external electromagnetic fields, was 
predicted as a direct evidence of the (not global but lo- 
cal) strong CP violation under very intense external mag- 
netic fields, and was observed in heavy ion collisions. Re- 
cently, an experimental evidence was presented by the 
STAR Collaboration at RHIC @. Since the discovery 
of the evidence, CME has been actively studied using 
non-perturbative techniques in QCD: P-NJL model [7], 
holographic QCD lattice QCD and so on. 

In this short note, in order to understand CME in 
QGP, we consider a generic linear sigma model (LSM) 
which is widely used as a key tool to understand the 
phase transitions. We find a universal mechanism for 
CME which is given rise to by a stable non-topological 
solitonic configuration of the (pseudo-)scalar mesons, so- 
called Q-ball [10( . Interior of the Q-ball is in the hadronic 
phase where the scalar mesons condense, while it is QGP 
outside the Q-ball. We find that the electric current along 
the external magnetic field arises in a similar manner dis- 
cussed in literature. In addition, as a consequence of the 
Q-ball, the electric current not only is a direct current 
but also has a small alternating current from quark mass 
terms. 

The generic LSM of the scalar mesons <&ij = q\q\^ is 
given by [l6[ 

C cS = Tr [<9 Al $9' I $ t - M(3> + - V^) 

+A(dct$ + det$ t ) , (1) 

where the metric is taken to be r)^ — diag(+, — , — , — ). 
The matrix M is proportional to the quark mass matrix 
M oc diag(m u , md, to s ), and the last term is a manifes- 
tation of the f (1)a anomaly in QCD. This Lagrangian 



enjoys the same symmetries as QCD. $ is singlet under 
SU{3) C x C/(1) B . The chiral SU(3) h x SU(3) R symmetry 
and U(1)a acts on <f> as <f> — !• e !Q [/ L $(/p, which are the 
exact symmetries if M = and A — 0. The currents 
corresponding to the axial part of these symmetries $ — > 
e ia U<PU, are given by J^ a = iTr[T a ($<9 M $t + M $t$ - 

fct^fc - where A a = 2T a (a = 1, . . . , 8) are the 

Gell-Mann matrices for flavor ST/ (3) and T° = -^=1.3. 

In order to discuss CME, we consider the electromag- 
netic field. The electromagnetic couplings of the quarks 
in QCD give rise to additional anomalies for the diago- 
nal elements of the axial currents, j5,a=o,3,8_ p ronl the 
effective theory point of view, these anomalies gener- 
ate couplings between the diagonal pseudo-scalars and 
the electromagnetic field. Our idea for the mechanism 
of CME is that non-trivial background for the diago- 
nal pseudo-scalars results in the electric current through 
these anomalous couplings. 

Motivated by such an idea, we concentrate on the di- 
agonal pseudo-scalars and the overall a field. Namely, we 
restrict $ as $ = diag($i, $ 2 , $3) = Ge ^' Ta ^ T ^^ T \ 
Then the Lagrangian £ e ff of the LSM, before coupled to 
the electromagnetism, is simplified to 

^(|a^i| 2 -^(*i+$;))-v(|#i| 2 )+A(n* i+ n*; 

i i i 

If Ma = A = 0, this model respects the U(l)'[_ R part of 
the axial J7(3)l-r symmetry. The corresponding current 
conservation laws d^ l J^ a = (a = 0,3,8) follow from 
the equation of motion, where the explicit forms of the 
currents are given by J^ a = 2iE i T?(<f> i <9 AI <f>* - <f>*9 A1 $i). 

Next we couple the electromagnetic field to this sys- 
tem. Since all the scalars we are considering are neutral, 
the only possible couplings are the anomalous ones ex- 
plained above. The explicit forms of such terms can be 
found by requiring that they should contribute to the 
anomalous current conservation law correctly. Our pro- 
posal is [TtJ 



U = 



^2 F ^ V ^(logSi-logS*) 
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HVpO p A T? 

c 1 fi- rL i'- L pa j 



(2) 
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where F^ v = \^ v P a F pa (e 



0t23 



I and each qi is the 



i-th element of the electric charge matrix for (u, d, s) 



2 



We have introduced P.. 



3-flavors, diag(|,— |,— 

(Po,P) = c^Im (J2i Qi log ®i) an d ignored the total 
derivative term in the right-most hand. With this term, 
the anomalous current conservation law is derived by us- 
ing the Euler-Lagrange equation for $i following from 
£ c ff + A, 



+V6iA(l[<Z> l -l[<S>*y a0 . 



(3) 



Once we restrict $ to the pseudo-scalar neutral mesons 
(77', 7To and 77), this reproduces the standard known form 
of the anomalous law. 

As mentioned above, the additional interaction £^ 
plays a role of an extraordinary source for the electro- 
magnetic field ((-1/2) A^J*)^'. 



3e 2 

47T 2 



(4) 



The Maxwell equations derived from the full Lagrangian 
including the Maxwell term, —jF^F^ + C e s + C-f, are 
modihed indeed, 



Vxi? 



dE 
~dt 



Jcm+J, V-E = Pem + J . (5) 



with Vx£+f=OandV-B = 0. Here J» = 

at em 

(Pcm, Jem) stands for the usual electric current, which 
vanishes for our model with the charge-neutral scalars. 
Thus the triangle anomaly in QCD results in the electro- 
magnetic currents. 

Note that these modified Maxwell equations are for- 
mally identical to those in the Maxwell-Chern-Simons 
theory [13] if the small fluctuation of 9 in [l|] is replaced 
by our LSM field $. Furthermore, relation between CME 
and axion strings and domain walls was studied in [l3[ . It 
has been proposed that CME occurs once the 9 parameter 
locally fluctuates 0, Q . Since current experiments sug- 
gest that the 9 parameter in the bare QCD Lagrangian is 
very small < 10~ 10 the origin of such a fluctuation is 
attributed to the effect of the medium in the QGP phase. 

In our approach with the generic LSM, on the other 
hand, CME is triggered by the Q-balls [l0| of the LSM 
field $, which are stable finite-size non-topological soli- 
tons. In the following, we shall show that the current of 
the CME is given by a typical frequency u attributed to 
the Q-ball, as will be found in (fTUj) . Note that our ar- 
gument is independent of the locally fluctuating 9 men- 
tioned above. 

In order to prevent inessential complexities, hereafter 
we will consider one-flavor model 

C = -\F^F^ + \d^\ 2 -V(\<P\ 2 ) + h{<S> + <i>*)-±A^J^ 

where J^ = — ^ s F v>1 d v \m (q 2 log3>) and h includes both 
the quark mass term M and the anomaly term A. This 



FIG. 1: A typical form of the scalar potential V. There is a true 
vacuum at a = (QGP). 
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FIG. 2: The short-dashed line is a typical form of U which allows 
the Q-ball. The solid line is the special case in which the two 
extremal values coincide, U(0) = U(<jq) = 0. 

Lagrangian has U(1)a symmetry if the last two terms 
vanish. Let us first construct the Q-ball in ft, = 
limit. We deal with the electromagnetic field as a back- 
ground field. [IH The Q-ball's charge, Q-charge, is the 
axial charge in this one-flavor model. 

The existence of Q-balls does not depend on the details 
of the system [10]. One requirement is that the scalar 
potential V(a 2 ) has a true vacuum at a — (QGP) as is 
given in Fig. [T] Let us make the following ansatz for a 
spherically symmetric Ty'-ball 

$ = cr(r)e l,? ' (t) , 



rf(t) = ut, 



(6) 



with r = The Euler-Lagrange equation for the profile 
function a leads 



ldU 
29^ 



0, 



U(p) = V(a 2 )-u 2 a 2 . (7) 



This system can be interpreted as a one dimensional clas- 
sical mechanics with the potential — U where r is "time" 
and a is "position." The term (2/r)a' plays a role of the 
damping force. Roughly speaking, the Q-ball is a soli- 
tonic solution connecting two extrema of —U(a). There- 
fore, the Q-ball exists when the extrimum at a ^ ap- 
pears and it is higher than that at a = 0, as is depicted in 
Fig. [5] The smoothness and finiteness of the solution re- 
quires a 1 = at both r — and r = 00. Coleman showed 
that a solution exists if cj is in the range Wq < uj 2 < fi 2 
(lOj . where fi specifies the curvature of the potential V 
at g = 0, = dV/da 2 \ a=0 , and w is Q-independent 
frequency, see below. 

The Q-ball can be best understood in the large Q- 
charge limit, where <j(r) resembles a smoothed-out step 
function. Then we assume that for small r less than a 
certain radius i?, a = const. > whereas for large r > R, 
a = 0, see Fig. [31 Namely, 



r<R 



r>R 



0. 



(8) 
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FIG. 3: A schematic picture of the spherical Q-ball. 

We ignore the contributions from the transition zone 
around r = R (surface of the Q-ball) which may be 
subdominant compared with the volume ones. In the 
limit, we obtain the energy E = lo 2 <j 2 V + V(a 2 )V and 
Q-charge Q — 2cjct 2 V with V being the volume of Q-ball, 
V = 47ri? 3 /3. A stable solution with fixed Q-charge is 
given by minimizing E with respect to three variables 
(a, uj,V) with the constraint Q fixed. First, E is ex- 
pressed in terms of (er, V) as E = \ + V{a 2 )V. Then 
by minimizing it with respect to V one gets E — Q-\J~^ 
and V = ; ^„ , which also determines cj in terms of a 

V4a 2 V 

as uj 2 (t 2 = V((J 2 ). Finally, we determine the value of a 
by minimizing E with respect to a. Let a ^ be the 
value of a for which E takes its minimum Eq, and luq be 
the corresponding frequency. Then, in summary, these 
values arc determined by 



Eo — Qw , 



= V'(a 2 ), 



«j%o%=V(o%), (9) 



where the prime stands for V'(a 2 ) = dV(a 2 )/d(a 2 ). The 
last two equations determine both <ro and ujq indepen- 
dently from Q. In fact these equations mean that the 
two extremal values of U(cr) = V(a 2 ) — uj 2 a 2 coincide 
and (To is one of the extrema with cr ^ 0, see Fig. [2] In 
this case, since er(r) can spend arbitrary long "time" r 
at the extremum a — oo, the solution can have an arbi- 
trary large volume V, and hence also the large Q-charge. 
Also, since the damping force in <[7|) is negligible after the 
long time, the profile o~(r) is well approximated by the 
smoothed-out step function. 

With the Q-ball at hand, we now see from Eq. (j4j that 
the electric current arises along a constant background 
magnetic field B (the electric field is assumed to be zero) 



J = 



3e 2 3 
UJB - 



(10) 



Here, we see that CME is a consequence of the existence 
of the Q-ball. The magnitude of our CME current is given 
dynamically by ui of the Q-ball. As found in ((9]), w is 
given by the LSM potential characterized typically by 
Aqcd- So it is natural to assume ujq ~ Aqcd- Using the 
expected value of the magnetic field eB ~ 10 4 [MeV 2 ] in 
heavy ion collisions [H , we obtain the magnitude of the 
CME current as | J\ ~ 10 5 [MeV 3 ] - 1(T 2 [fur 3 ]. 



Also, natural size of the Q-ball is ~ Q 1/3 [fm], if all the 
dimensionful parameters are of order Aqcd ~ l[fm -1 ]. 

Note that this important frequency w plays a role of the 
so-called chiral chemical potential (A5 [5(. The chemical 
potential can be introduced through the change do$ — > 
(do + i/J,5)Q in (J6j> . The relevant terms for this change 
arc the ones with time derivatives: 



Ido^l 2 - d Im(log$ 



3e 2 q 2 



A ■ B 



8tt 2 

(5 Im(log$) + Ai 5 ) 



-A ■ B. (11) 



8tt 2 

Taking the Q-ball solution $ = a(r)e lut in the theory 
without fj>5 is equivalent to considering a static solution 
$ = u(r) in the theory with /Z5, if we identify ^5 = lj. 
In the latter case, the current (ITU1) is supplied by the last 
fi 5 A- B term in ([TT]). 

Since the Q-charge is preserved, this non-topological 
soliton is fairly stable. However, the h($ + $*) term in 
Eq. ^ breaking explicitly U(1)a may destabilize the Q- 
ball, which would result in destroying a constant supply 
of the electric current. So, let us next analyze the effect 
of h($ + <&*). We expect that, if h is sufficiently small, 
U(1)a is broken only weakly and Q-ball still lives long. 
In the following, we shall derive a condition for h to have 
the stability, and find that Q-ball are fairly stable, but 
with a new interesting feature of alternating CME cur- 
rent component. 

We treat h as a small parameter and we expand fields 
with respect to a small dimensionless parameter e as 

h 



1 ' 

do Wo 



(12) 



where $0 is the Q-ball solution in e — > limit. Again, 
we consider the large Q-ball limit given in Eq. (0). We 
would like to solve the equations of motion 

F = d^<S> + $y' - h = 0, (13) 

order by order in e as F = F + eFi + e 2 F 2 + ■ ■ ■ =0. 
The zeroth order Fq = is the Q-ball equation which we 
have solved. As we have seen, this gives us ujq = V'(o-q). 
The next-to-leading order is F\ = with Fi given by 

5^$ 1 + $ir(|$ | 2 ) + $o($o*i + *o$D^"(l*o| 2 )-^o 
We solve this equation in two region, r < R and r > R, 
separately. In the former region we put 4>o = croe iw °* 
while $0 = in the latter region. Therefore, the next- 
to-leading order solution is given by 



\r<R 



3ai — 1 
3ai 



1 



7^0 



3ai 



rfoe 



2iu)ot 



$0 



aQto 2 /V'(Q), with a\ 
-leadi: 

-36 (2 + a 2 



v '("a) 



next-to-next-to-leading order is readiliy solved by 

-l' 



(14) 



The 



a 2 



9ai 



\r<R 



9(2- 



8 (3ai + 2) 
-6ai + 2 



-foe 



-iujQt 



2(3ai 



-a e 



iujot 



3(-3ai + 2) 
'8(3ai + 2) : " 



;o- e 



3iu)Qt 
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and $9 



r>R 



0, where we have defined ui 



The mass of the Q-ball up to this order is evaluated as 



E = Quj + - 



3 (-3(5 + 2a 2 " 1 )af + 3ai + 2) 
(3ai + 2) 2 



Qw e 2 



where we have shifted the origin of the energy in such a 
way that the energy density outside the Q-ball becomes 
zero, and Q = 2uo<JqV as before. Note that the first 
contribution starts at the order of e 2 and the next is of 
order e 4 . Therefore, the variation in energy is negligible 
if e is sufficiently small, namely h <C uqoJq. By using 
WqCT 2 = ^(ctq) ( see ©), this condition can be rewritten 
as hao <C V(o"q). This condition is quite natural since 
it just means that the perturbation h($ + $*) is small 
compared to the original potential I^d^l 2 ). Hence we can 
expect that the Q-ball is stable against the perturbation. 

Interestingly, a contribution of order e arises in the 
electric current as 



J 



3e 2 



q'Bojo 1 - e 



3ai-2 
: 3ai+2 



coswo^ ' 



(15) 



Thus the quark mass term and the U(1)a anomaly in 
QCD eventually give rise to a small alternating CME 
current. This is a new feature of CME by the Q-ball. 

Let us finally make a comment on a possibility of the 
Q-balls by the other pseudo-scalar mesons such as ir° and 
■q. In addition to C/(1)a, there exists an axial part of the 
chiral symmetry S{7(3) L _ R e S*C/(3) L x SU(3) R in QCD. 
It is straightforward to construct the Q-balls by using a 



t/(l) 2 _ R subgroup in SU(3)l-r, for instance, a pionic 
Q-ball with the T 3 generator. Since the i7(l) 2 _ R is also 
anomalous by the electromagnetic interaction, the CME 
current arises as in the case of the ?/-ball [20| • 

In summary, we present a useful formalism of LSM 
which can explain CME in QGP via the non-topological 
soliton, Q-balls. The electric current arises along the 
external magnetic field and it has a small alternating 
current as a consequence of quark mass terms and the 
£/(1)a anomaly in addition. Furthermore, the interior of 
the Q-ball is the hadronic phase. So we predict that there 
may be a lump of hadrons, Q-ball, in QGP. It might have 
a certain contribution in the cooling process of QGP and 
hadronization. 

In cosmology, the Q-ball is thought of as a candidate 
of so-called boson stars [15j . We hope that our study 
may open a new direction to create baby boson stars at 
RHIC, LHC and FAIR. 
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